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Abstract We briefly summarize properties of quantum Hall states with a pairing 
or clustering property. Their study employs a fundamental connection 
with parafermionic Conformal Field Theories. We report on closed form 
expressions for the many-body wave functions and on multiplicities of 
the fundamental quasi-hole excitations. 1 
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1. Introduction: new states of matter in the 
quantum Hall and BEC arenas 

The fractional quantum Hall effect has unveiled states of matter that 
can be characterized as incompressible quantum fluids with topological 
order. Such states are formed in a two-dimensional electron gas, at very 
low temperature and in the presence of a strong perpendicular magnetic 
field. It has been recognized early on [17, 7] that the excitations over 
fractional quantum Hall states obey fractional braid statistics: a con- 
figuration of N quasi-holes over a fractional quantum Hall ground state 
forms a one-dimensional representation of the braid group Bn, where 
the braiding of two quasi-holes is typically represented by e ia7T , with a a 
rational but non-integer number. The requirement that particle states 
have to represent the braid group rather than the permutation group is 
special for two dimensions: the braid group is the fundamental group of 
the configuration space of identical particles only in two dimensions. On 
general grounds it is known that, for two-dimensional quantum systems, 
higher dimensional representations of the braid group are allowed 
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(see [13] for an early reference). In such a situation, the braiding of par- 
ticles is represented by matrices, and since matrices in general do not 
commute, this leads to the notion of non-abelian statistics. 

It is now believed that the 'non-abelian statistics scenario' is realized 
in novel types of quantum Hall states, which are characterized by a 
pairing or clustering of electrons under quantum Hall conditions. There 
exists concrete experimental [32] and numerical [21] evidence that the 
simplest of these states, the 'pfaffian state' proposed by Moore and Read 
[22], exists in nature, as the state of matter underlying the quantum 
Hall effect at filling fraction v = |. At present, this example stands as 
the lone confirmed 'sighting' of a non-abelian quantum Hall state. One 
expects that paired and clustered quantum Hall states can be realized 
in multi-layer quantum Hall systems with sufficiently strong interlayer 
tunneling, or in situations where the spin of the electrons is not fully 
polarized. 

It has recently become clear that quantum Hall states with clustering 
correlations are relevant for a very different class of physical systems, 
namely Bose Einstein condensates (BEC) of cold atoms in a rotating 
trap. In a regime where the rotation frequency Q is of the order of the 
frequency uj set by the confining potential, and where the healing length £ 
is very large, a rotating BEC becomes formally analogous to a quantum 
Hall system of bosonic particles. In numerical studies, the quantum 
liquid that is formed in a rotating BEC at vorticity v = Afb oson /iV vortex = 
k/2, has been identified with the Read-Rezayi quantum Hall state with 
order- k clustering [10]. 

In this paper we briefly review the properties of a variety of paired 
and clustered quantum Hall states, with particular attention for states 
formed by spinful electrons. We present closed form expressions for 
ground state wave functions and discuss the multiplicities and statistics 
properties of quasi-hole excitations. Throughout our presentation, we 
stress the role of parafermionic conformal field theory as a crucial tool 
in the theoretical description. 

2. The fractional quantum Hall effect 

The discovery of the fractional quantum Hall (fqH) effect [28] was 
truly remarkable and unanticipated. At fractional filling fraction v = 1/3 
a quantum Hall (qH) plateau was observed. The filling fraction is defined 
as the ratio of the number of electrons and the number of available states 
in the lowest Landau level: v = N/N^, where is the number of flux 
quanta piercing the sample, and iV the number of electrons. Soon after 
the discovery, Laughlin made a fundamental break-through by proposing 
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his by now famous wave functions, which describe the qH effect at filling 
fraction v = ^, where m is an odd integer [20] 

*z\z 1 ,..., ZN ) = i[(z l -z 1 r . (i) 

i<j 

Here and below we display reduced qH wave functions ^f(z), which are re- 

~ |- |2 

lated to the actual wave functions ^(z) via ^f(z) = ^(z) exp (— J2i -$r) 

with I = \pj^ the magnetic length. 

Although the qH effect occurs at relatively high magnetic fields, it 
was soon realized that the electron spin can indeed play an important 
role. The spin-polarized Laughlin states were generalized by Halperin, 
who proposed a set of spin-singlet wave functions [18] 

$™ +1 ' m+1 ' m (4,...,4;z},...,4)= (2) 

11(4 - 4) m+1 11(4 - 4) m+1 11(4 - 4) m > 

i<j i<j i,j 

where zj and zj are the coordinates of the spin-up and spin-down elec- 
trons, respectively. The state Eq. (2) has filling fraction v = 2/(2m + l). 

2.1 The qH effect-CFT connection 

Bulk connection. Following Moore and Read [22] one observes 
that it is natural to view (lowest Landau level) qH wave functions as 
conformal blocks of electron-type operators in a suitable chiral conformal 
field theory (CFT) in 2+0 dimensions. This point of view is related 
to the fundamental role of Chern-Simons field theories for qH systems 
(compare with [33], where an explicit link between Chern-Simons theory 
and CFT is established). 

As an example, the Laughlin ground state wave function (1) is ob- 
tained as 

*l = lim {^) mN2 {V e { Zl ) . . . V e (z N ) : e-^^oo) :) , (3) 

Zoo ^OO 

with V e (z) =: exp(iy / my}) : a chiral vertex operator in the c = 1 chiral 
CFT describing a single scalar field ip compactified on a radius R 2 = m. 

Edge connection. While bulk excitations over a qH fluid are gapped, 
one expects gapless excitations at the edge of a sample. Following Wen 
[29], one observes that the edge excitations are described by a chiral 
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Luttinger Liquid or chiral CFT in 1+1 dimensions. In the example of the 
v = ^ Laughlin states, one again has the scalar field theory at R 2 = m. 
The neutral operator p = iy/md(p is associated with edge density waves, 
while vertex operators of type V q (z) =: eKp{iq^/mip) : represent charged 
edge excitations, the charge being equal to qe with — e the charge of the 
electron and ^ the charge of the fundamental quasi-holes. 

2.2 Fractional statistics in the fqH effect 

In the case of an abelian qH state, changing the magnetic field by 
one flux quantum $ = | results in the creation of a quasi-hole (or 
particle, depending on the sign of the change). These quasi-holes can 
have remarkable properties, such as a fractional charge. Also, the quasi- 
holes over the Laughlin fqH states are anyons, i.e. they realize fractional 
braid statistics [17, 7]. The fundamental phase for the braiding of two 
such excitations is given by e ! ». 

Closely related to this are the fractional exclusion statistics of these 
same excitations [16, 19]. Focusing on edge excitations, one can show 
that the gapless, charged edge excitations of an abelian qH state satisfy 
a form of exclusion statistics closely related to that of Haldane [16]. A 
particularly natural choice of basis for the edge excitations employs edge 
electrons and quasi-holes as the fundamental quanta [11]. In this basis, 
the exclusion statistics parameter matrix is diagonal with self-exclusion 
parameters equal to m (for the edge electrons) and ^ (for the edge 
quasi-holes) . 

For general abelian qH states, one may argue [1] that the statistics 
matrix G of edge excitations (in a specific basis) is of the form 

G = K e ©K e ~ 1 , (4) 

where K e is the so-called i^-matrix that characterizes the topological 
order of the qH state (see for instance [30]). 

3. Paired and clustered qH states 

Prompted by the observation of a qH effect at filling fraction v = | 
[32] a number of novel qH states have been proposed. Among these 
is the Moore-Read (MR) state or pfaffian state, which is characterized 
by a p-wave pairing of the electrons [22, 14]. A generalization, where 
the pairing is replaced by a clustering of order k was proposed by Read 
and Rezayi (RR) [25]. In [5] two of the present authors made a further 
generalization to a class of spin-singlet qH states also characterized by 
a clustering into fc-plets of electrons. 
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3.1 Spin-polarized states 

The wave function of the (spin polarized) MR state is given by 

*MR(*i) = Pf(— ) IK* " ^') M+1 , (5) 
Zl Z 3 i<j 

where Pf (M i:j ) = 2N/ 2( N/2)l E CT sgn(a) Ur=i M *(2r-i),<r(2r) is the pfaffian 
of an antisymmetric matrix. For the wave function to be antisymmetric 
(we are describing electrons) M needs to be an odd integer, which implies 
an even-denominator filling fraction v = jjti ■ 

The MR wave function can be viewed as a correlator in a c = 3/2 
CFT, consisting of a free scalar field and a Majorana fermion. The 
electron operator becomes ip(z) : exp(iy/M + 1 tp c ) : (z). The correlator 
of N electron operators (and a suitable background charge) splits into 
a product of vertex operators, giving the Laughlin part of the wave 
function, and a product of fermion fields, which gives the pfaffian factor. 

Upon generalizing the Majorana fermion to the Zfc parafermions [34] 

associated to the coset s "ffi fc , one obtains the clustered states of [251. 

u(l) L J 

Their wave functions are constructed in the same way as the MR wave 
function, with explicit parafermion factors brought in by the electron 
operator. The result is a state in which the electrons form clusters of 
order k rather than pairs. The filling fraction takes the form v = fcA ^ , 2 , 
with M an odd integer. Note that for k = 1 the Laughlin states (with 
m = M + 2) are recovered, while k = 2 gives the MR states. 

In ref. [9], Cappelli et al. presented a particularly simple construction 
for fc-clustered qH states, and they showed that the resulting wave func- 
tions are equivalent to the ones originally obtained by Read and Rezayi 
[25]. In this construction, one starts with an abelian state of k types of 
discernible electrons, with the same total filling fraction as the clustered 
states. For the RR states, this abelian state takes the form [9] 

^b=nnte a -4) 2 > ( 6 ) 

a=l i<j 

where a labels the different types of electrons. In the construction, we 
will assume that M = 0; the wave function for M / are obtained by 
multiplying the final result with the Laughlin factor IIj<j( z j — z j) M ■ Note 
that the wave function (6) can be thought of as k copies of the clustered 
state with parameters k = 1,M = (which is a bosonic Laughlin state 
at filling fraction v = \). The clustered states are obtained by effectively 
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making the electrons indiscernible, which is achieved by symmetrizing 
the expression (6) over all the electrons. The wave function for the RR 
states can now be written in the following form 



niK 

a=l i<j 



~a\2 



i<j 



(7) 



where it is assumed that after the symmetrization procedure all the 
electrons are labeled by a single index. 

The clustering property of the wave function (7) for M = is inherited 
from the abelian state (6). The clustering property of order k means 
that the wave function does not vanish if up to k particles are brought 
to the same position; however, if k + 1 particles are at the same position, 
the wave function does vanish. To see this, note that k electrons all 
of different type are brought to the same location, the wave function 
(6) does not vanish. The symmetrization makes sure that this property 
holds for any k electrons of the clustered state (7) (with M = 0). 

3.2 Spin-singlet states 

The clustered states discussed in the previous section do not include 
a spin degree of freedom. As was shown in [5] and [2], the clustering 
property can be extended to spin-singlet states of spinful electrons, lead- 
ing to what have been called 'non-abelian spin singlet' (NASS) states. 
As stated before, the Halperin states are spin-singlet analogues of the 
Laughlin states. In the same way, the NASS states of [5] are spin-singlet 
analogues of the RR states. To construct their wave functions, two bo- 
son fields are needed (for charge and spin) in addition to the 'higher 

rank' parafermions associated to j^jyp [12]. The NASS states have fill- 
ing fraction v = 2 kM+3 w ^ n M an odd integer. For k = 1 the Halperin 
states (2), with m = M + 1, are recovered. Explicit expressions for the 
NASS wave functions of [5] can be found in [4]. The form of the wave 
functions presented there is similar to the explicit RR wave functions in 
[25]. Alternatively, a formulation similar to the construction for the RR 
states presented above is also possible. We again start with k copies of 
the k = 1, M = states, which are, in this case, the Halperin states with 
labels (2,2,1). Again, this is a bosonic state. To obtain the fermionic 
states, one has to multiply in the end with an overall Laughlin factor 
Ui<j( z i ~ z j) M w ith M odd. The abelian state from which the clustered 
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spin-singlet states can be obtained takes the form 

* ab = n ( U(4' a - 4* a ) 2 IT(4' a - 4 ,a ) 2 Il(4' a - 4 a ) ] • (8) 



Again, we have k types of electrons, but now also the spin is taken into 
account. It is not so hard to find the filling fraction of this abelian 
state, which is given by v = which is indeed the filling fraction of 
the NASS states with M = 0. To obtain the clustered state, one has 
to symmetrize the wave function (8). However, because we are dealing 
with spin-singlet states, the symmetrization procedure is more involved 
than for the RR states. The conditions which are to be satisfied by 
the spatial part of the wave function in order be a spin-singlet are the 
Fock cyclic condition, and the conditions that the wave wave function 
be symmetric in the spin-up and spin-down coordinates separately. In 
general, this can be achieved by the following (generalized) symmetriza- 
tion procedure (see [4] for more details) . First antisymmetrize in z\ , z\ , 
followed by antisymmetrization in z^z^, etc. After the antisymmetriza- 
tion of zjj^ 2 , z^/ 2 i the wave function has to be symmetrized over all the 
spin-up coordinates, and finally, over the spin-down coordinates. This 
gives rise to clustered spin-singlet wave functions. Though we do not 
have a proof, we believe that via this procedure, we obtain wave func- 
tions which are equivalent to the ones described in [4]. For the case at 
hand, the antisymmetrization steps as described above do not seem to 
be necessary (as was the case for the explicit formulation in [4]). Thus 
we arrive at the following form for the wave functions of the clustered 
spin-singlet states 



*NASS 












k 


IM ,a 


-4* a ) 2 lM' a 






a=i \ 




i<j 




x Y[( Xi - x_ 


i) j 






i<j 











(9) 



t i 

where X{ stands for either z\ or z\. <S 2 t ;2 j. denotes symmetrization over 
the spin-up and spin-down particles. The clustering property of these 
spin-singlet wave functions can again easily be read off from this formu- 
lation. For more details, see [4]. 

Recently, another type of clustered spin-singlet qH states was pro- 
posed [2]. These states are based on parafermions corresponding to 
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another rank 2 affine Lie algebra, namely the j^yp parafermions. Also, 
two chiral boson fields for the charge and spin degrees of freedom are 
present. At level k = 1, the only parafermion field present in this theory 
is the Majorana fermion, which also appeared in the theory for the MR 
state. It is thus to be expected that the state corresponding to the level 
k = 1 affine Lie algebra sb(5) is related to the MR state. This is indeed 
the case; the wave function also contains a pfaffian factor [2] 

*&, = Pf(-^-)^ +1 ' M+1 ' M (4^j) > (io) 

where \j>^+ 1 > M + 1 ' M j s th e Halperin wave function (2). For M odd, 
this describes a qH system at filling u = 2 m+i ■ This wave function 
is different from the previous ones in a few respects. First of all, the 
(fundamental) excitations over this state show a separation of the spin 
and charge degrees of freedom (SCsep). This is a consequence of the 
structure of the underlying Lie algebra symmetry, as was pointed out 
in [2]. In effect, there are spinon like excitations, with spin-i and no 
charge, and holons which carry charge 2 m+i without spin. 

The other feature which is different in comparison to the clustered 
states (7) and (8) is the clustering property. The clustering of the 'spin 
charge separated state' (10) (with M = 0), is in fact a clustering of the 
spin-up and spin-down particles separately. Note that the wave function 
for M = has poles when spin-up particles are at the same location 
as spin-down particles. Nevertheless, we will discuss the clustering for 
the case M = 0. In the case k = 1, up to two particles of the same 
spin can be brought to the same location while the wave function is still 
non-zero. If we first bring z\,z\ — > z^ and z\,z\ — ► z^, and then send 
z^ — > the wave function remains non-zero. The clustering is thus 
different from the clustering of the NASS states (8), which vanish when 
any k + 1 particles are brought together. 

We refer to a forthcoming paper [3] for further analysis of the spin- 
charge separated qH states for general k. 

4. Quasi-holes over paired and clustered qH 
states 

In a BCS superconductor, where electrons are paired up, the funda- 
mental flux quantum is reduced to = A. The same phenomenon 
occurs in the paired and clustered qH states, and this means that in- 
serting a single flux quantum $o creates more than a single quasi-hole. 
For the spin-polarized states of [24] the number of quasi-holes is given 
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by n = kANtf), where AN^ is the number of excess flux quanta . For the 
spin-singlet states of [5], this relation becomes m) + = 2k AN^, were 
n}'^ denotes the number of spin-up and down quasi-holes, respectively. 

The quasi-holes over the paired and clustered qH states carry frac- 
tional charge and satisfy non-abelian braid statistics. They can be stud- 
ied with the help of the associated CFT. The wave functions of states 
in the presence of quasi-holes are obtained by inserting into the CFT 
correlators the appropriate quasi-hole operators, which consist of a ver- 
tex operator part and a spin field of the parafermion theory. In the 
case of the MR state, this is the spin field a of the Ising model and the 
quasi-hole operator becomes a(w) : exp( 2 /J^ +1 fc) ■ (w). 

The non-abelian statistics have their origin in the non-trivial fusion 
rules of the parafermion spin fields. In general, there is more than one 
way to fuse the fields in the correlator to the identity; for n spin fields 
the number of ways will be denoted by d n . The braiding of n quasi-holes 
is then represented by a matrix of size d n x d n . 

4.1 Braid statistics 

The simplest example that exhibits the non-abelian braiding is the 
situation where 4 quasi-holes are added to the MR state. In this case 
di = 2, so there are two distinct states which, following [23], we write 
as $(4qh,o) and ^(4qh,i)_ g tarting f rom t h e state ^(4qh,o) ; and braiding 

two of the particles, we find the following transformation 



Wave functions for the MR state with n quasi-holes can be written as 



^(4qh,0) _^ £.1 ^(4qh,0) + ^(4qh,i) 




(11) 



[24] 



*MR,qh(zi, ■ ■ ■ ,Z N ;Wl,-- -,W n ) 



1 



Y[(zi - Zj) 



M+l 



2 (N-F)/2( N _ F y 2 \ 



i<j 



F {N-F)/2 



®{Zg(F+2l-l),Zg(F+2iy, Wl, ...,W n 
(z a (F+2l-l) - Z a (F+2l)) 



) 



X 



E^IKS) n 



o-eSjv k=l 1=1 



(12) 
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where 



$(z 1 ,z 2 ;w 1 , . . . ,w n ) = 



((n/2)!)2 



1 



n/2 

~ W r(2r~l))(z2 ~ W T (2r)) ■ 



rSS„ r=l 



(13) 



The integers mi, . . . , must be chosen such that they satisfy < mi < 

rri2 < ■ ■ ■ < mp < — 1 , giving rise to a degeneracy = ^ J, ^ • The 

number F is interpreted as the number of unpaired electrons in the 
excited state. 

The braid matrices for n quasi-hole excitations were obtained by 
Nayak and Wilczek [23], who showed a direct connection with the ro- 
tation matrices of the group SO(2n). We refer to [27] for more general 
results on braid matrices. 

4.2 Quasi-hole counting formulas 

The CFT approach to the excited state wave functions and their braid 
properties is highly efficient. One would like however, to 'keep both 
feet on the ground' and understand the fundamental degeneracies that 
characterize the non-abelian statistics in a more direct way. This can be 
done by selecting a (ultra-local) hamiltonian that has the qH state as its 
ground state, and then (numerically) studying the spectrum of excited 
states. 

These numerical computations are most easily performed by studying 
a small number of particles in a spherical geometry. By tuning the 
number of flux quanta to the value = — S, where S is the so- 
called shift [30], one realizes the qH state as the unique ground state. 
Cranking up the number iV^ and performing a numerical diagonalization, 
one obtains characteristic degeneracies for quasi-hole excitations. 

Following [24], we first explain the counting of degeneracies for the 
case of the MR state. To understand the degeneracies of quasi-hole 
excitations on the sphere, two effects should be taken into account. The 
first is a choice of fusion path or, equivalently, a choice of numbers 
F and mi, . . . ,mp in the formula (11). The second effect is the so- 
called orbital degeneracy: the quasi-holes are not localized on the sphere, 
but can occupy one of a finite number of available orbitals, each of 
which is characterized by a definite angular momentum L z . These orbital 
degeneracies are well-known from the analysis of integer and abelian qH 
states. 

For the MR state, the orbital degeneracy factor depends on the num- 

( F*\ 

ber F of unpaired electrons. Fixing this number F, we have d n different 
choices for the quasi-hole wave function. To each of those we can asso- 
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/ H=E + n \ 

ciate [24] an orbital degeneracy factor equal to I 2 I . Putting 

it all together, we have the following total degeneracy 

W«)=E(f )(^ n +n ) . (14) 

in agreement with numerical results [24]. The degeneracies d n relevant 
for a situation where n quasi-holes are at fixed positions are recovered 
by suppressing the orbital factors, d n = X^drf^ = 2 n / 2_1 , where the 
sum is over even (odd) F for N even (odd). This number is in agreement 
with a direct count of the number of fusion paths of n Ising spin fields 
[23]. 

For the more general clustered qH states, the degeneracies have ba- 
sically the same form: an orbital part and an intrinsic part, stemming 
from the non-trivial fusion rules. The difference is however, that we can 
not rely on explicit wave functions to handle the intrinsic degeneracy. 
One can work around this by extracting from the parafermion CFT the 
relevant combinatorial factors, using the methods put forward in [26, 8]. 

For the RR states, the counting has been worked out in [15], with the 
result 

1 / N—F , x 



#(AT,n;fc)=E{;}^^V "J , (15) 

with n the number of quasi-holes, n = kAN^. The symbols { F }k rep- 
resent the degeneracies due to the fusion rules. In [8, 15], these were 
described in terms of recursion relations; explicit formulas (based on bi- 
nomials) for general k can be found in [6]. The sum d n = which 
equals the total number of fusion paths for the spin fields contained in n 
quasi-hole operators, sets the dimension of the braid matrices for braid- 
ing 2 out of the n quasi-particles. 

For the NASS states of [5], the counting goes along the same lines, 
with the additional complication that we have to deal with two spin 
components, which are combined in a non-trivial way. This is reflected 
in the counting formulas by a doubling of the number of binomial factors. 
By inserting an amount AN^ of extra flux, one creates quasi-holes, which 
can have either spin. The total number of quasi-holes is fixed, w) + = 
2kAN ( p. The symbols {}k now depend on four parameters { F p 2 }k and 

we have d n ^ n i = J2f 1 ,f 2 {f 1 F 2 }k- The case k = 2 has been worked out 
in detail in [4], where the results are checked against numerical data. 
Explicit results for the symbols {}& can be found in [6]. The counting of 
quasi-hole degeneracies over the spin-charge separated spin-singlet states 
will be discussed in [3]. 
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The numbers d n (for both spin-polarized and spin-singlet states) are 
easily extracted from the known fusion rules of the su(2)k and su(3)k 
CFTs. For both su(2)s and su{3)2 the numbers d n are Fibonacci num- 
bers. The asymptotic behavior for n — > oo is found to be 

d n ~[2cos-^]" (16) 

for the RR states, and 

2-7T 

4,„[l + 2co8— ]* (17) 
for the NASS states, where p = w) + . 

4.3 Exclusion statistics and K matrix structure 

In [1] a proposal was made for a K-matrix structure of the paired and 
clustered qH states discussed in this paper. It was established that the 
exclusion statistics of edge excitations over these states (in a suitable 
basis) can be captured by a statistics matrix of the form (4), supple- 
mented by a prescription that some of the particles described by this 
matrix be viewed as pseudo-particles. We refer to the first paper of [1] 
for a physical picture underlying these ET-matrices, and to the second 
paper of [1] for mathematical details. 
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Notes 

1. Note that we adopt a slightly different notation than the one used in [24, 15]. Here, n 
denotes the number of quasi-holes, rather than the number of excess flux quanta. 
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